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Abstract 

M. Porto (Phys. Rev. E 63 (2001) 030102) suggested a system consisting of Coulomb interacting 
particles, forming a linear track and a rotor, and working as a molecular motor. Newton equations 
with damping for the rotor coordinate on the track x, with a prescribed time-dependence of the ro- 
tor angle 0, indicated unidirectional motion of the rotor. Here, for the same system, the treatment 
was generalized to nonzero temperatures by including stochastic forces and treating both x and 
via two coupled Langevin equations. Numerical results are reported for stochastic homogeneous 
distributions of impact events and Gaussian distributions of stochastic forces acting on both the 
variables. For specific values of parameters involved, the unidirectional motion of the rotor along 
the track is confirmed, but with a mechanism that is not necessarily the same as that one by Porto. 
In an additional weak homogeneous potential field U(x) =const-x acting against the motion, the 
unidirectional motion persists. Then the rotor accumulates potential energy at the cost of thermal 
stochastic forces from the bath. 

PACS numbers: 05.30.-d,05.70.-a 
Keywords: Porto system, molecular motor 



I. INTRODUCTION 



There are many constructions of so called ratchet systems devised to turn stochastic 
perturbations or noise acting on classical or quantum systems into a unidirectional motion. 
Perhaps the most famous example is the 'wind- mill' system analyzed by Feynman De- 
tailed analysis mostly revealed that such systems could work but in experimental set-ups 
devised to violate the Second law of thermodynamics, they fail jl], |2j . This fully corresponds 
to results of multiple theoretical analyses. 

Surprisingly, the uni-directionality of the motion has recently been found to be determined 
not only by a type of asymmetry of the (say) ratchet potentials incorporated but rather 
dynamically 0, |4], The activity in this direction resulted into construction of a relatively 
realistic model consisting of a molecular rotor on a linear track, both consisting of particles 
interacting by just Coulomb forces ||. Under a prescribed time-dependence of the rotor 
angle 0(t) entering the potential energy of the system, the classical Newton equation for its 
linear coordinate x along the track yields that the rotor unidirectionally moves along the 
track. The prescribed time- dependence 0(t) may also be stochastic what introduces the idea 
that one can easily convert the system, upon adding a potential U (x) slightly increasing in 
the direction of the rotor motion, into a machine converting energy of the stochastic influence 
of the bath, that might be of thermal origin, into a mechanical potential one. This would 
mean so called perpetuum mobile of the second kind and a violation of the Second law of 
thermodynamics. 

Porto explicitly refuses this possibility as far as his calculations were concerned - see 
a note in this respect on page 3 of 0. Our more general calculation below, however, 
confirms the above conjecture. Here, we should like to add that this fact is by no means 
objection against the Porto system and its analysis using the classical Newton equations. 
Instead, our treatment rather generalizes that one of 0. First, one can easily justify adding 
stochastic forces to the Newton equations, converting them into the Langevin ones. In this 
respect, our calculations below provide a proper extension of those by Porto || to finite 
temperatures. The Langevin equations have well understood quantum counterparts in the 
Mori or Tokuyama-Mori identities f?|, |8|, [| what makes the calculations even more relevant. 
Second, the rotor angle B has not been here ascribed any specific time-dependence like 
in || but has been assumed to change according to another Langevin equation with the 
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same temperature entering the stochastic force correlations functions; this is rather more 
physical. In fact, one should not be very much surprised by the challenge to the Second law 
provided by this amended treatment of the Porto model generalized to finite temperatures. 
The reason is twofold: 

• First, there is at present increasing evidence, both theoretical [[H], [□], [T^, [L3|, |T^, 



15 1 and experimental |T6| , fL7j , [18j as well as combined one [HI, in favor of potential 



violability of the Second law. 

Second, the present model and its present treatment, though they are perhaps more 
realistic than often in similar situations, are still just theoretical and their relation to 
Nature may be not as obvious as it might seem at the first sight. 



II. FORMULATION OF THE PROBLEM 



In contrast to the Porto study, we impose no prescribed time-dependence of the G(t) 
variable but describe both the linear coordinate x(t) and the G(t) angle as two general- 
ized coordinates describing the rotor (its position as well as the rotation) with their time- 
dependence determined by a pair of the corresponding coupled Langevin equations. This 
method supersedes and definitely surmounts the Newton equations used by Porto, in at least 
the sense that 

• the 0(t) variable has its time-dependence determined also from a dynamic equation, 
and 

• this approach respects existing connections between dissipation (friction) incorporated 
and by its effect decisive in the Porto model, and properties of the stochastic forces 
on the right hand side of the Langevin equations that were completely ignored in ||. 

So, the basic set of equations reads 

mx + V x = -|-$(:r, 9) + T x (t), JQ + k& = ~$(x, 6) + T e (t). (1) 
ox a<~) 

Here m, J, rj and k are the rotor mass, its moment of inertia, and the linear and angle 
friction constants. As for the stochastic forces T x (t) and re(t), they should fulfil standard 
equations 

(T x (t)T x (t')) = 2 v k B T K 5(t-t'), (r e (t)r e (f )) = 2 K k B T K 5(t-t'), (T x (t)T e (t')) = (2) 
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ensuring, in absence of the potential $>(x, 0), that the asymptotic mean squared velocities 
((x) 2 ) and ((G)) 2 ) fulfil the equipartition theorem with temperature Tk- Here (. . .) designates 
the ensemble average. The stochastic forces are here represented as 

r„(t)=i; r e (t)=£/f5(t-tf), (3) 

i i 

where tf and tf are statistically independent times of impact events. Times between two 
such succeeding impacts are exponentially distributed 

w x {t) = i exp(-t/Q, w @ (t) = i exp (-*/*©) (4) 

t x lQ 

(t x and te being the mean waiting times) while distributions (f) and w^ & (f) of (also 
stochastic and statistically independent) impact 'forces' ff and ff are, on grounds of phys- 
ical arguments, assumed Gaussian. For @ to be satisfied, these distributions should read 



V A7rr)t x k B T K y ^Kt e k B T K 

In order to get the model fully defined, on must specify the potential energy <3>(x, 6). It 
is connected with the definition of the system as given by Porto |J. The track consists of 
charges q > at tW = {(0.1 + n)6, 0, -0.256} and -g at tH = {(0.4 + n)6, 0, -0.256} with 
n being arbitrary integer. So, the track is neutral, periodic but without inversion symmetry. 
The rotor has the total mass m and apart from potentially other neutral atoms, it consists 
of four point charges q^ v = (— 1) 1+A V, /i, v = 1 or 2. Designating as above position of the 
rotor center of mass as x, their positions are = {x, 0, 0} + r^sh^G + AQ^ V ), 0, cos(0 + 
AG^)}, AG^ = (v - h/2)tt. The two radii are n = [0.005 + 0.045Max(- 7 , 0)] and 
r 2 = [0.005 + 0.045Max(7, 0)]. The parameter 7, — 1 < 7 < +1, determines the 'gear' of the 
Porto molecular motor model. With that, the potential energy $(x, O) of the configuration 
of the charges reads 



2 +00 
$(x,9) = gg' £ (-1) 1+ " £ 

/j,i/=l n=— cx 



1 1 



t { n +) -x^(x,9)| |tl } -x^(a;,e)| 



(6) 



Calculation of the potential $(x,0) requires some care because the individual sums 
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III. NUMERICAL TREATMENT 



Before numerical treatment, it is useful to reformulate the problem in dimensional units. 
We introduce dimensionless coordinates and time by 



Then ([I]) reads 



where 



x{t) = X{T) -b, t = t x ■ T. 



^x(T ) + E±x { T) = -^(x,e) + lx{n 
|,e(T) + ^e(^-.e>,e) + 7 e(T) 



c _ V ,- r _ K - _ qq'(t x ) 2 _ qq'(t x ) 
& — — z xi — —Ixi •> X — — > •> 8 — 7T~ 

m J mb 6 Jo 



(7) 



(8) 



{i,v=\ n=— o< 

Formulae (^|) are then replaced by 



|x(+) 



x»„(xb,e)\ |d- } -x^(x6,e)| 



7x(T) = £i^(T-lf), le{T)=Y,F?*{ T -Ti 



(9) 



(10) 



where, instead of (||), we have distributions of the (statistically independent) dimensionless 
times T„ x and T® 



W X (T) = exp(-T), W e (T) = ^exp(-T^). 



(11) 



As for the distributions of the (again statistically independent) 'forces' F x and F®, we have 

W pX {F) - v ,,-,., <'xp 



—F 



2 (5m 2 b 2 

HQ 3 



W pe {F) 



(3J 2 



exp 



—F 



f3J 2 



, (3 = \/{k B T K ). 



(12) 



4rei e (4) 2 . 

With that, the problem has now been solved numerically, using a random number generator 
for simulating the stochasticity involved. The input data used were 

£ = 10, K = 0.001, T x = 13.79, JF = 0.1092, 



0.0101, 



(3J 



127.5 



4£(Q 2 '4/Ci/e '"" (13) 

what could correspond to, e.g., m = 100 x the proton mass, b = 10 -6 cm, t x = 10~ 10 sec, 

£e = 10~ 8 sec, temperature Tk = 300 .K" and J = 10mb 2 . 
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IV. NUMERICAL SOLUTION 



Series for ip(X, 0) in (|9]) as well as those for its partial derivatives converge very slowly. 
For numerical solution, however, it is necessary to calculate as accurately as possible. Using 
extrapolation methods, we were able to calculate all these functions reasonably fast with pre- 
cision of 12-13 digits at least. Rotor-track potential energy (p(X, G) is a periodic function in 
both variables, with period 1 in X (b in x) and tt in Q. Form of this function for the 'forward 
gear' 7 = +1 is shown in Fig. 1. Form of the potential profile for 7 = —1 ('reverse gear') 
can be deduced from Fig. 1 using identity (p(X, 0,7 = —1) = —(p(X, + 7r/2,7 = +1). The 
respective positions of extremes for 7 = +1 and their values are X = .1179500143638, = 
.175365440363, (p = -.33522791779058053399124 (minimum) and X = .3802499856363, 
9 = -.175365440363, cp = .33522791779058053399124 (maximum). For 7 = -1, the values 
are X = .3802499856363, = 1.395430886432, up = -.33522791779058053399124 (min- 
imum) and X = .1179500143638, = 1.746161767158, up = .33522791779058053399124 
(maximum; all the digits being valid). 

Lengths of waiting times between two succeeding impact events were generated from uni- 
formly distributed random numbers r and r' G (0, 1) using relations W X (T) = — ln(r) and 
W 6 (T) = — £§-ln(r'). The Gaussian distributed impact forces F t x and F® were generated 



using procedure gasdev from [20 



Between any two impact events, set of equations (||) was solved using the Bulirsch- 
Stoer method (procedure bsstep from with numerical constants as in ([13]). Values of 

constants Jj^y and 4K ^ t - p are then .0101 and 127.5, respectively. The ratio |j = 0.01. 
Initial values of X and were always chosen at minimum points of the potential energy ip 
(see above); initial values of ^ and ^ were always chosen as zero. 

Because of the nonlinearity of (§), one must expect a chaotic behavior (extreme sen- 
sitivity to numerical values of initial conditions as well as finite accuracy of the numerical 
calculations). Fig.'s 2a and 2b illustrate the onset of the chaos. For our accuracy 10~ 13 in in- 
tegration of the differential equations (§) and the above accuracy of calculation of derivatives 
of uj, the chaos starts to appear always at about T = 1200. 



6 



V. RESULTS 



Fig.'s 3 and 4 illustrate time-dependence of solution for X(T) for different sequences of 
the random numbers (i.e. impact events as well as the values of the impact forces) involved. 
In agreement with Porto's results, the tendency of behavior of X(T) is to increase (for the 
'forward gear') and to decrease (for the 'reverse gear') with increasing T. One of the curves 
in Fig. 3 in interval (1,1000) provides an exception from the rule. Nevertheless, extending 
the calculations up to T = 4000, general tendency of final increase was anyhow finally found 
to appear also here. Because of lack of the random forces in the Porto calculations, no such 
exceptional behavior could be found by him. 

Fig. 5 shows a new remarkable feature of our calculations. As Porto, in ||, has in 
fact solved the problem for just X with a prescribed dependence of 0(T), and because 
his calculations involved no stochasticity, his only mechanism for the forward and backward 
motion of the rotor along the track was that one owing to the tendency of X(T) to relax fast, 
at any time, to a minimum of the potential <p(X, O) with the prescribed instantaneous value 
of 0. Thus, in his case, coordinate X(T) continually moves along valleys of the potential 
profile, being driven by the prescribed time dependence of 0. Because we have in our case 
also the stochastic impacts on the coordinate X included, the rotor motion along the track 
could also be owing to the impact-induced hops of the rotor between two neighboring minima 
of tp, without any essential change of the 'slow' variable 0. (In the whole interval (0,25) in 
Fig. 5, changes by just about 1/27 of a turn.) Fortunately, as Fig. 5 shows, direction of 
the prevailing motion owing to this new mechanism is the same as that one owing to the 
©-induced motion of X(T) by Porto. Thus, the two mechanisms of the rotor motion along 
the track rather support each other. 

The whole beauty of the behavior of the system just reported becomes explicit when we 
realize that the unidirectional tendency of the rotor motion could survive even when we add 
a smooth potential acting against the above two mechanisms. In order to make the the 
situation more clear, we have added, to the potential <f(X, 0) above, a linear potential 

8<p(X,e) = c-X, c = 0.01. (14) 

In Fig. 6, we can see some results of the simulation. The curves correspond to otherwise 
exactly the same parameters as in Fig. 3. Worth noticing is that even the curve, that 
virtually corresponded to the backwards motion of the rotor in the limited time-interval 
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investigated, turns now up to describe, like all other curves, motion of the rotor against the 
homogeneous potential field QI^). Thus, because of the periodicity of ip(X,Q), the system 
starts now to accumulate the potential energy S(f(X, 0). The question is what is the source 
of this acquired potential energy that could be later used in any prescribed way. We can 
see its only source - this is the above stochastic behavior simulating influence of a real 
bath on the system. If so, the behavior of course becomes incompatible with the standard 
thermodynamics. This illustrates how the physics behind, in particular, the Second law of 
thermodynamics is still little understood. 

VI. CONCLUSIONS 

The generalized Porto model introduced above, as defined by the Porto system || and the 
Langevin equations ([I]) describing its behavior in external stochastic fields modelling influ- 
ence of the thermal bath, provides a relatively realistic model that makes numerical studies 
of its time development reasonable. Though the description by the Langevin equations is 
classical, it should be taken seriously because of existing quantum counterparts of these 
equations. This makes the basic conclusions obtained relevant. These are, in particular: 

• Confirmation, at finite temperatures, that the rotor has a tendency to move along 
the track in prevailingly one direction only. This tendency was already found, in his 
simplified treatment corresponding inter alia to zero temperature, by Porto ||. 

• Identification of another mechanism (in addition to that one found by Porto himself) 
also leading the unidirectional motion of the rotor along the train. 

• Survival of the behavior even when the motion goes against a weak potential field. 

The latter characteristic behavior indicates possible incompatibility of the model and its 
description via the Langevin equations with principles of the statistical thermodynamics. 
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VIII. FIGURE CAPTIONS 



Fig. 1: Plot of the dimensionless rotor-track potential energy ip for 7 = +1 as a function 
of X and 6. 

Fig. 2a: Time dependence of angle 6 as a function of T for 7 = +1 and the same values of 
all the input parameters as specified in the main text. In all three cases, the impact 
times as well as the random forces were identical. Only in time intervals between any 
two individual impacts, different integration time steps (1/32, 1/16, 1/10, 1/8 and 
1/7 for curves a) to e), respectively) were used. The onset of chaos appears at about 
T = 1200. At T = 1000, the values of 9 obtained still coincided to 5 digits. 

Fig. 2b: The same for X as a function of T for the same values of all the input parameters. 
At T = 1000, the values of X obtained still coincided to 6 digits. Notice that all the 
curves increase with increasing time. 

Fig. 3: Dependence of X on T for five different sequences of the random numbers involved 
for the 'forward gear' 7 = +1. 

Fig. 4: The same for the 'reverse gear' 7 = — 1. The same sequences of random numbers 
as in Fig. 3 were used. 

Fig. 5: Details of time-dependence of X(T) (7 = +1) for curve 'c' of Fig. 3 before the first 
impact event in 6 occurs. Times T of the first six impacts in X were 4.4506, 5.0489, 
5.9216, 9.7589, 9.9581, and 10.7401. The first impact in 6 appears at T = 216.8751. 
Worth noticing is how the impacts cause, because of the periodicity of <p, sudden 
changes of X = x/b (on the vertical axis) by integers. 
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